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MODELING OF PLASMA ETCHING IN MICROELECTRONICS 



V. P. Derkach, I. P. Bagrii, 
and G. A. Chechko 



Computer simulation and experimental results are used to establish the possibility of efficient control of 
plasma etching rate and nonuniformities. 

Plasma etching is a basic technology for the manufacturing of integrated circuit topologies. The available models 
[1-6] do not fully reflect the macrokinetics of the plasma etching process and are not suitable for modern equipment. 
The problem of correct modeling of the etching process therefore remains highly relevant and important in many appli- 
cations. 

The following processes take place in a plasma-chemical reactor: production and recombination of chemically 
active particles (CAPs) in the volume, convective and diffusive transport of CAPs, chemical reactions between CAPs 
and the etched material, heterogeneous recombination of CAPs. In a general form, the plasma etching process can be 
described by a system of differential equations 

9 — gradp + n&v + f, (1) 

div r= 0, (2) 

■~ + (wgrad) c, = 0,Ac, ± Q t „ (3) 

where p is the gas density, v is the flow velocity vector, p is the pressure, p. is the viscosity of the gas, f is the volume 
force, C; is the concentration of CAPs of species i, D; is the diffusion coefficient of CAPs of species i, Q oi is the capacity 
of the volume source responsible for production or recombination of CAPs of species i. 

Translated from Kibernetika, No. 5, pp. 27-33, September-October, 1990. Original article submitted December 
21, 1988. 
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Fig. 1. Schematic diagram of the individual treatment 
reactor: 1) perforated anode, 2) plate, 3) cathode. 
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Fig. 2. Distribution of CAP concentration (a) and etching 
rate (b) over the plate radius for k = 2.4 and various h 0 : 1 ) 
h 0 = 0.1, 2) h 0 = 0.2, 3) h 0 = 0.3, 4) h 0 = 0.4, 5) h 0 = 0.5, 
6) h 0 = 0.6, 7) h 0 = 0.7, 8) h 0 = 0.8, 9) h 0 = 0.9. 



The Navier— Stokes equation (1) describes the hydrodynamic motion of the gas in the reactor. The continuity 
equation (2) is an expression of the law of conservation of matter. Equation (3) describes convective diffusion during 
plasma etching. Note that the parameter Q oi describes only volume processes. Particles produced or destroyed on the 
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Fig. 3. Distribution of etching rate over the plate radius 
for various k and h 0 [a) k = 2.4, h 0 = 0.5, b) k = 2.4, h 0 = 
4.1, c) k = 0.19, h rt = 0.5, d) k = 0.03, h 0 = 0.2] and vari- 
ous values of the parameter A: 1) A = 0, 2) A = 0.1, 3) 
A = 0.2, 4) A = 0.3, 5) A = 0.4, 6) A = 0.5, 7) A = 0.8. 
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surface are allowed for by boundary conditions. The boundary conditions for the system of equations (l)-(3) have the 
following form: 

v = 0 on all surfaces, (4) 
Ddcj/dn = kcj m on all surfaces, (5) 
Cj is finite at the reactor center, (6) 

where n is the direction of the normal to the surface, k is the rate constant, of the surface reaction (chemical reaction 
or recombination process), m is the kinetic order of the reaction. In physical terms, expression (5) describing the surface 
reaction kinetics implies that the flow of CAPs on the surface is determined by surface reactivity. 

Solving the system of equations (l)-(6) is a fairly complex undertaking. It can be simplified if we exploit some 
specific features of the plasma etching process. Calculations show that the Peclet number, which characterizes the 
efficiency of diffusion and convective processes, is 0.05-0.3 [1}. For comparison, the Peclet number of a liquid is 10 4 . 
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Fig. 4. Schematic diagram of the reactor in the PLAZMA 
1251 facility: 1) quartz cone, 2) HF inductor, 3) plasma 
production zone, 4) perforated screen, 5) treatment zone, 
6) plate, 7) cathode, 8) body. 
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Fig. 5. Etching nonunif ormity of polysilicon (1-5), silicon 
nitride (6), and silicon oxide (7) in the PLAZMA 1251 
facility in a SF 6 + 20% 0 2 mixture at 160 Pa pressure, 40 
W HF power, for various perforated anode designs from 
Table 1: 1) anode No. 5, 2) anode No. 4, 3) anode No. 3, 
4) anode No. 1, 5) anode No. 2, 6, 7) anode No. 5. 



Plasma etching therefore may be regarded as a quasifixed system. Moreover, since the etching time is much greater than 
the CAP lifetime, plasma etching also may be viewed as a quasistationary system. We also assume that CAP production 
and recombination is uniform throughout the volume, so that the system of equations (l)-(6) reduces to 

DfACi — k Pl c\ -\- G t = 0 (7) 

with the boundary conditions 

DjdCj/dn = kjCj m on all surfaces, (8) 

C; is finite at the reactor center, (9) 
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where k p .C; { and G l are the homogeneous recombination rate and the production rate of CAPs of species i, k p . is the 
homogeneous recombination rate constant, t is the kinetic order of the reaction. In general, the coefficients D v k p ., kj 
may be variable. 

Equation (7) with boundary conditions (8)-(9) is the basic macrokinetic equation of plasma etching. Solving this 
equation, we obtain the distribution of the CAP concentration, determine the derivative dc jdn, and calculate the etching 
rate v = ^ j ; . = £ £> . (dc,/dn) in flow units [cm- 2 sec -1 ]. In the reactor this equation describes the loading eff ect, 

on the plate it describes the rate, selectivity, and nonuniformity of etching, and on an element of the circuit topology 
it describes the etching profile. 

Let us consider the limits of applicability of the proposed model. The diffusion approximation is valid when the 
CAP free path A is much less than the characteristic dimensions of the system. For plasma etching in the pressure range 
10-100 Pa, A « 100-10 Mm. For the reactor and the plate, these conditions are satisfied; for a topology element, A ~ d. 
Equations (7)-(9) therefore must be used with care when describing plasma etching processes on a topology element. 
If A > d, then the diffusion approximation is inapplicable to the description of plasma etching on a topology element 
and we have to consider free-molecular motion of the CAPs. 

Solving equations (7)-(9) for a specific plasma-etching reactor by numerical methods, we obtain ail the process 
characteristics, but this is a complicated and time-consuming process. Specializing the plasma-etching conditions and 
the boundary conditions of the problem, we can apply the general model to obtain a number of interesting practical 
results. As an example, consider the plasma etching of silicon-containing materials (Si*, Si 3 N 4 , Si0 2 ) in an individual 
treatment installation with SF 6 + 0 2 plasma. 

The recent increase of plate diameters has encouraged the transition to individual treatment facilities. In some 
facilities, the plasma discharge zone is separated by a metallic grid from the processing zone. The schematic form of 
the plasma reactor in such systems is shown in Fig. 1. For definiteness, we refer to the lower electrode as the cathode 
and to the upper electrode as the anode. The etching plate is located on the cathode; the anode is perforated. This anode 
design makes it possible to investigate etching nonuniformities for various concentration profiles of the CAPs arriving 
at the plate surface. 

Since CAPs are generated only in the plasma discharge zone, etching in SF 6 + 0 2 plasma is performed by CAPs 
of one species only — fluorine atoms. The chemical etching reactions are first-order reactions in CAP concentration and 
CAP recombination occurs only on the etching surface [8]. Assuming that the CAP diffusion coefficient is independent 
of CAP concentration, we obtain from the general equation (7) Ac = 0. 

Initially, let us consider the determination of etching nonuniformity in the absence of CAP sources on the anode 
surface. The diffusion equation in cylindrical coordinates for the axially symmetric case has the form 
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with the boundary conditions 



■ -gj-C, U) = A. c (/-,<)), (11) 

where c 0 is the given concentration at the reactor inlet (it is determined experimentally), h is the interelectrode distance, 
r„ is the plate radius. 

We introduce a new concentration function C(r, z) = c — c 0 . Problem (10), (11) is transformed to the form 



C(r 0 ,z)=0, <j£.(r.h) = 0, 

D~ (r,0) = kC(r,Q)+kc 0 . (13) 

The solution of problem (12), (13) is best sought in dimensionless form: 

d 2 C | 1 dC i d 2 C __ q 

d? r dr dz- ' (14) 

C(l,z) = 0, (r, /i 0 ) = 0, (r, 0) = 

dz dz (15) 

= ~kC (r, 0) + k, 

where C = C/c 0 , f = r/r 0 , z = z/r 0 , k = r 0 k/D, h 0 = h/r 0 , 0<r<l,0<z<h. The general solution of the problem 
obtained by separation of variables is representable in series form 

C (r, z) = V [A n ■ Sh Mn (h 0 - ~z) + fl n -Shu„i] (!-t n r"), (16) 

where A n = n n 2 and p n (n = l,2,...,oo) are the roots of the equation J 0 (/i) = 0. The series coefficients A n and B n are given 
by the expressions 



.4,, = 8 n -ch\jL,Ji n , (17) 
B n = 2k/{ixl ■ J 1 ((!„)•[ 1 — ch ji.A, (ch n„/t 0 + (18) 
+ kShn n h 0 )]}, 

where J 0 (/i), Jj(^) are Bessel functions. 

Let us write dimensionless expressions for the CAP concentration at the anode surface for z = h 0 and for the flow 
on the etched plate for z = 0, which is equivalent to the etching rate: 



C (r, h 0 ) = £ B n Sh a n h a -J 0 {\i n -r) + 1, (19) 
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(20) 



We used formulas (19) and (20) in some computations, which revealed a marked decrease of CAP concentration 
in the direction of the OZ axis (Fig. 2a). Indeed, as the CAPs progress by diffusion from the edge of the plate to its 
center, the CAPs are gradually absorbed by the chemical reaction on the etching surface. 

It follows from (19) that if a gas mixture is additionally forced through the anode surface, supplying a 
supplementary concentration distribution density near the anode according to the formula 



the CAP concentration for z « h„ will be maintained equal to c 0 . This, however, may be insufficient for ensuring 
uniform etching with acceptable tolerance over the entire surface of the plate in the full range of parameter variation 
in practice. Such uniformity apparently can be ensured only if the CAP concentration at the anode surface is distributed 
by the law 



where the function g(f) is unknown and, by axial symmetry of the problem, satisfies the conditions g'(0) = 0 and g(l) - 
0. Thus, uniform etching can be achieved only if the additional inflow of gas mixture through the perforated anode 
modifies the plasma discharge so that the following supplementary CAP concentration is established on the anode at z = 
h 0 : 



Calculations show that the function g(r) must be positive. 

In general, the problem of improving the etching uniformity and finding the function F(r) can be solved 
mathematically, by specifying a constant CAP flow at the plate surface for z = 0 (constant etching rate over the plate). 
In dimensionless variables, the problem takes the form 



f 6 = -^B„.ShnA-4(w), 



C(?,h) = c a + gCr), 



Fg (r) = g (r) - J B n -Sh » n h 0 .JJvJ). 



d' 2 c , 1 dc t d 2 c 
~dr r ^TUF^ ~W = ' 

c(r 0 ,z) = c 0 , c(r,ft) = /(r), D ~- (r, 0) =q = kc lt 



(22) 



(21) 



where c l is a constant. 

We solve problem (21), (22) by separation of variables. The solution is represented by the series 




(23) 
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where A n and B n are given by the expressions 



B n = -j— r fr./WX 

^•/?Cu n )Sh^/ij U J 

X J 0 (^ry r -c a rl-J 1 (|i B )/|* B j, 

r§4(n n )Sh 2 (-^-ft Lo J 
X [2*. j dr - c^J, (u,)/u.„J 2 ^ ir ° 



We need an additional condition in order to find the unknown function f(r). Isolate a cylinder with cross-section 
radius r 0 within the active part of the reactor and define the CAP flow through all the surfaces of this cylinder: Q x is 
the CAP flow through the cylinder wall, Q 2 is the CAP flow through the perforated anode surface, and Q 3 is the CAP 
flow to the etching plate. By the law of mass conservation, we clearly have Qj + Q 2 = Q 3 and after simple manipulations 
we obtain in expanded form 

For a given flow q = kc x , the unknown function f(r) can be found numerically from the integral equation (24). When 
the flow q, and therefore the etching rate, are variable, the function f(r) will also vary. Numerical solution of Eq. (24) 
involves considerable difficulties and requires a special mathematical apparatus, because the problem is ill-posed. This 
means that there may exist several functions f(r), all of which satisfy Eq. (24), and the function f(r) determined by the 
numerical method is not necessarily the unique solution that corresponds to the physical content of the problem. 

The problem of improving the etching uniformity is thus better solved by an approximate fitting of the function 
f(r). From physical considerations we conclude that it has a nearly parabolic form: f(r) = c 0 + A(r 0 2 — r 2 ) = c 0 + g(r), 
where A is an unknown constant. Varying A, we can obviously ensure approximate constancy of the flow (etching rate) 
on the plate surface. 

The diffusion equation and the boundary conditions in this case have the form 



c (r 0> z) = c 0> c (r, h) = c a + A {r* — r% 
D~(r, 0) = kc (r, 0). 



Introducing a new concentration function C(r, z) = c(r,z) — c 0 and changing to dimensionless variables, we transform 
problem (25), (26) to 

1 dC &C 

^ 2 + r d ~ r ^ «fe 2 ' (27) 
C(U) =-0, C(r,h B ) = A(\-?), 
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(vg) 



where A = Ar 0 2 /c 0 , z e [0, h 0 ], r £ [0, 1], k > 0. 

Problem (27), (28) is solved by separation of variables and its solution is represented by series (16) with the 
coefficients A n and B„ given by the expressions 



Calculation results regarding etching nonuniformity for various parameters h 0 , A, and k are presented in Fig. 
3. Analysis of the curves indicates that etching nonuniformity can be altered by forcing an additional gas mixture into 
the interelectrode space and a CAP distribution density can be chosen for which uniform etching is achieved on the 
entire plate surface with practically acceptable tolerances. The graphs in Fig. 3 show that as the interelectrode distance 
increases, the supplementary CAP flow rate forced through the perforated anode must be increased for plasma-etching 
control. Calculations show that as the coefficient k decreases, the etching uniformity improves, which also facilitates 
the process of control. 

Theoretical modeling of etching rate and nonuniformity leads to a number of important practical conclusions. 
First, a uniform distribution of CAP concentration on the anode for etching of high-k materials does not ensure uniform 
etching. Second, by shaping the CAP flow on the etching surface, e.g., by means of a perforated screen, we can control 
the etching uniformity between wide limits. 

Experimental tests of the theoretical results were made on a PLAZMA 1251 individual treatment facility. Films 
of polysilicon (0.55 /im thick), silicon nitride (0.15 |im), and silicon oxide (0.1 fim) on silicon plates 100 mm in diameter 
were etched in a SF 6 + 20% 0 2 mixture at pressures of 160 Pa and 40 W high-frequency power. The reactor is shown 
schematically in Fig. 4. The plate treatment zone is separated from the plasma production zone by a perforated metallic 
screen which, on the one hand, protects the plate from bombardment by plasma ions and, on the other, shapes the CAP 
flow hitting the surface of the plate. The chamber design is very close to the theoretical reactor shown in Fig. i and is 
ideally suited for studies of nonuniform etching under individual treatment conditions: cathode diameter 120 mm, 
reactor height h 0 - 30 mm, plate diameter 100 mm (r 0 = 50 mm). 

Chamber design and etching technology are developed so as to ensure a uniform concentration of CAPs on the 
screen surface, and the screen (the anode in our notation) is usually perforated uniformly with constant transparency 
along the radius. In practice, however, this chamber design does not produce a uniform distribution of CAP con- 
centration in the treatment zone. The solutions obtained above may be used to produce both uniform and nonuniform 
CAP concentrations. By varying the anode transparency according to a certain law, we can produce the desired CAP 
distribution profile and with it the desired etching uniformity. 

The following method was used to obtain varying anode transparency. The anode radius was partitioned into 
three zones. Each zone was characterized by a certain grid spacing and hole diameter. By varying the zone radius, the 
grid spacing in the zone, and the hole diameter, we could regulate the magnitude and the distribution of anode 
transparency between wide limits. The anode characteristics used in our experiments are listed in Table 1. 

Comparison of the theoretical model with experimental results should allow for the fact that the calculations were 
carried out for certain values of the dimensionless variable k, which is a function both of the discharge and of the 
interaction of the material with the CAP — the rate constant k of the chemical reaction. The exact value of the rate 
constant was not available in most cases, although it is well known that k si > k SisN4 > k Si02 [9]. Therefore, only a 
qualitative comparison of the theoretical and experimental values was possible. Polysilicon films have a high rate 
constant for the reaction with fluorine atoms, the main CAP in the SF 6 + 0 2 discharge. The experiments with these films 
therefore produced the most significant results. 



A u = ( Ch a n h 0 + kSh \i n h o r l -{B n — 2*/(i4-y i (n n )]}, 



{2A/[^ n -J 1 {ii n )Sh^ n h 0 ]} x 




Etching of polysilicon films with a standard anode of uniform transparency (1 in Table 1) produces a higher 
etching rate at the center of the plate than at the periphery (Fig. 5, curve 4). This etching pattern is determined by 
nonuniformity of the initial CAP flow reaching the anode surface in the chamber of our design (Fig. 4). The observed 
etching pattern could not be changed by varying the technological conditions (pressure, mixture composition, HF 
power). Application of a perforated anode with the same distribution of transparency (3 in Table 1) enhances the etching 
nonuniformity. Anodes in which the distribution of transparency is the inverse of the CAP distribution produced by 
the chamber design (2, 4, 5 in Table 1) alter the etching pattern all the way to complete inversion: the etching rate at 
the periphery of the plate is higher than at the center (Fig. 5). Similar results, although less sharply pronounced, were 
also obtained for silicon nitride and silicon oxide etching (Fig. 5). 

Experiments have shown, however, that the use of perforated anodes in systems of this kind is in itself 
insufficient to ensure etching uniformity of polysilicon films which meets modern technological requirements. The 
existing reactor design does not produce a sufficiently uniform CAP flow along the chamber diameter. It is also 
necessary to ensure a uniform CAP flow on the anode surface by altering the reactor design. 

Note that the numerical values of the plasma etching kinetic constants (CAP diffusion coefficients, etching rate 
constants, heterogeneous and homogeneous recombination rate constants, CAP production rates) and their dependence 
on technological process parameters have been little studied [6]. As a result, our model ignores these factors and, while 
correctly estimating the qualitative picture, it does not pretend to produce an exact description of the process. Further 
development of the model must account for a number of additional factors, such as the effect of plasma discharge on 
CAP concentration and transport during etching, ambipolar diffusion of high-energy particles, convective transport of 
CAPs and high-energy particles, electron and ion bombardment, UV radiation. 

Plasma etching, being a heterogeneous chemical reaction which includes CAP production, recombination, 
transport processes, and chemical reactions on the surface, may be described to first approximation by the CAP diffusion 
equation of the form DjAc; — kp.C;' + G; = 0 with the boundary conditions of the form D ; grad c ; = kft™ reflecting the 
kinetics of surface reaction. The proposed model describes heterogeneous chemical reactions (on the assumption that 
the kinetic constants k i? k p ., G;, D i5 1, m, and other factors dependent on external plasma discharge parameters are 
known) and reflects the macrokinetic features of the plasma etching process. The kinetic constants entering the model 
may be obtained from experimental data or from special theoretical studies. In application to individual treatment of 
plates, the model parameters are the geometrical dimensions of the reactor and the plate, the rate constant of the 
chemical reaction of the CAPs with the material, and plasma discharge parameters (diffusion coefficient and con- 
centration of CAPs). 

Our modeling has shown that in order to ensure uniform etching of materials with high chemical rate constants, 
a nonuniform CAP flow must be produced on the etching surface. A suitable concentration pattern can be created by 
a perforated anode of variable transparency. The modeling results are supported by experimental etching data for 
polysilicon, silicon oxide, and silicon nitride in the PLAZMA 1251 individual treatment facility with SF 6 + 20% 0 2 
mixture at 160 Pa pressure and 40 W high-frequency power. 
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METHODS TO FIND THE MOST PREFERRED ELEMENT ON A SET OF 



Stochastic interactive procedures are considered, based on paired comparison of efficient solutions with 
parametrization of the Pareto-optimal set. The convergence conditions of these procedures are analyzed. 



Multicriterion problems often arise in connection with mathematical modeling of economic, technical, social, 
and other complex system. Let us consider one of such problems. 

Assume that the set D of all feasible solutions is given in E n . Each (not necessarily feasible) solution x is an m- 
dimensional vector of criterion values f(x) = (fi(x),...,f m (x)), m>2, which is improved by increasing any of the criteria, 
while keeping the other criteria fixed. It is required to find a solution x* € D which corresponds to optimal (in a sense 
to be defined later) criterion values. This problem will be written in the form 



Given problem (1), a solution x which for some y £ D satisfies f;(x) < f ; (y) and f k (x) < f k (y) for at least one 
1 < k < m is a priori not the optimal solution x*. Eliminating all the solutions x from the feasible set D, we obtain the 
set D of Pareto-optimal (efficient) solutions. If D is closed and bounded and the functions f s (x) are continuous, then the 
set D is clearly nonempty and at least one of the solutions of the problem max xeD fj(x) belongs to this set. 

In various applications, the set D usually consists of more than a single solution [1]. In some problems of the form 
(1) we have [2] D = D, and additional assumptions are needed in order to determine x*. In what follows we assume that 
a decision maker exists who is capable of comparing any two solutions (in particular, any two solutions from D) by the 
corresponding criterion values and choosing that which is better. The search of the solution x* of problem (1) is an 
interactive procedure by which the decision maker chooses the best solution from D. The decision maker compares the 
criterion vectors generated by computer; the function of the computer in this case is to make sure that these vectors 
represent feasible (and, moreover, efficient) solutions. 

The mathematical model of the decision maker in this interactive procedure is a binary relation R defined on 
the set of all m-dimensional vectors f as follows: f(*)Rf( 2 ) if the decision maker judges the solution with the criterion 
value vector fl 1 ) not worse than the solution with the criterion vector f( 2 ). Problem (1) now can be restated as the 
problem of finding a solution x* € D which for any x e D satisfies f(x*)Rf(x). The purpose of this paper is to describe 
methods for solving this problem under appropriate assumptions about the properties of D and R. 

Let D be a convex compactum in E n and let the relation R be complete, reflexive, transitive, and continuous [3]. 
The distinctive feature of the methods considered in this paper is that all the solutions presented to the decision maker 
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